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基于 改进 的 NSGA2 算法 考虑 病 患 公 平 性 及 


医院 运作 成 本 的 病床 配置 优化 研究 
陈 ， 龙 ， 刘 勤 明 '， 叶 春明 ， 李 佳 翔 


(上 海 理 工大 学 管理 学 院 ， 上海 200093) 


摘 要 : 针对 目前 医院 病床 调度 存在 运营 成 本 较 大 以 及 医 患 关系 之 间 公 平 性 的 问题 ， 提 出 一 个 考虑 医院 运作 成 本 和 

病 串 公平 性 下 单 科 室 病 床 分 配 的 多 目标 随机 规划 模型 。 首 先 ， 根据 基于 医院 的 相关 政策 ， 提出 一 个 考虑 响应 性 与 准 
入 性 的 权重 测度 指标 来 反映 医 患 关系 的 公平 性 ， 并 考虑 医院 的 运作 成 本 建立 多 目标 随机 规划 模型 ; 其次， 为 方便 算 

法 求解 ， 引 入 线性 化 方法 将 复杂 模型 处 理 成 混合 整数 线性 模型 ; 最 后 ， 采 用 改进 后 的 NSGA2 算法 对 多 目标 问题 求 

解 ， 并 对 算 例 进行 不 同 的 数值 实验 ， 通过 调整 不 同 的 参数 进行 相对 应 的 灵敏 度 分 析 ， 改进 后 的 算法 提升 了 算法 的 收 

伊 性 与 多 样 性 ， 实 验 结果 验证 了 模型 的 有 效 性 和 适用 性 
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Bed allocation based on improved NSGA2 considering patient equity and hospital costs 


Chen Long, Liu Qinming', Ye Chunming, Li Jiaxiang 
(Business School, Shanghai University of Technology, Shanghai 200093,China) 


Abstract: To address the problems of large operating costs and fairness between doctors and patients in the current hospital 
bed scheduling, this paper proposes a multi-objective stochastic programming model for single department bed allocation 
considering hospital operating costs and patient fairness. First, this paper proposes a weighted measure of responsiveness and 
access to reflect the fairness of the doctor-patient relationship based on relevant hospital-based policies, and establishes a 
multi-objective stochastic programming model considering the operational cost of the hospital; second, to facilitate the 
algorithm solution, a linearization method is introduced to process the complex model into a mixed integer linear model; 
finally, the improved NSGA2 algorithm solve the multi-objective problem The improved algorithm enhances the 
convergence and diversity of the algorithm, and the experimental results verify the validity and applicability of the model by 
adjusting different parameters for different numerical experiments. 


Key words: improved NSGA2; hospital wards; equity; multi-objective optimization problem; equity measure index 


0 ”引言 采用 蒙特 卡 洛 法 仿真 求解 。 如 Song 等 人 根据 病 患 在 医院 中 
的 分 布 情况 ， 设 计 了 基于 排队 论 的 急诊 室 就 诊 系统 ， 系 统 会 
在 医院 诊疗 阶段 ， 病 人 需 依据 自身 病情 接受 诊疗 ,然而 根据 病人 的 诊疗 结果 进行 病床 分 配 ， 有 效 地 绥 解 了 
病人 在 无 序 下 地 排队 通常 会 导致 医疗 资源 的 浪费 。 其 中 ， 病 堵 的 情况 外 。 周 雄伟 等 人 为 了 解决 患者 在 就 诊 过 程 中 等 
床 资源 作为 医疗 资源 中 的 关键 环节 ， 能 否 合理 配置 直接 影响 。 时 间 过 长 ， 排 队 拥堵 等 方面 的 问题 ， ee 
了 医院 的 收容 水 平 、 运 作 水 平 、 服 务 水 平等 。 考 虑 到 资源 有 门诊 资源 调度 的 最 优 策略 四。 上 述 的 研究 主要 集中 在 医院 单 
限 性 ， 内 部 可 使 用 面积 以 及 运营 成 本 等 问题 ， 如 何 对 病床 资 ”一 科室 的 病床 分 配 问 题 ， 然 而 不 同 科室 的 病人 就 诊 情况 差异 
源 优化 成 为 一 个 亚 需 解决 的 问题 赔 。 医 院 作为 国家 医疗 卫生 ” 较 大 ， 通 常会 导致 菜 一 科室 的 病床 闲置 ， 而 男 一 科室 的 病床 
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体系 的 基础 ,还 需 考虑 不 | Ee dei 医院 。 紧张 的 情况 。 因 此 在 提升 医院 病床 利用 率 方面 ， 单 一 对 某 一 

在 日 常 运作 过 程 也 需要 控制 成 本 t。 因 此 ,在 保证 患者 就 诊 。 科室 的 病床 分 配 往 往 效率 不 高 。 

公平 性 的 同时 也 要 考虑 到 人 去 作成 本 情况 成 了 难点 。 随 着 学 者 们 对 病床 资源 分 配 的 进一步 深入 研究 ， 多 目标 
针对 医院 的 病床 分 配 优化 问题 ， 国 内 外 学 者 们 提出 了 许 ”优化 方法 也 成 为 学 者 们 求解 多 科室 分 配 优化 问题 常用 的 方法 

多 不 同 的 解决 方法 。 目 前 病床 分 配 问题 主要 分 为 单一 科室 的 之 一 ， 但 是 医疗 资源 的 复杂 性 经 常 导 致 多 个 目标 函数 之 间 产 

病床 资源 优化 四 以 及 多 科室 或 考虑 其 他 因素 的 联合 优化 。 对 ” 生 冲 突 0% 习 ， 对 此 国内 外 学 者 根据 病床 资源 的 特性 提出 了 许 


于 前 者 ， 诸 多 学 者 主要 采用 排队 论 以 及 其 衍生 的 相关 模型 来 。 多 目标 分 配 优化 策略 。GreenM31 等 指出 医院 在 病床 分 配 的 过 
解决 诸如 手术 室 、 重 症 监护 室 、 术 后 观察 室 等 单 科室 病床 的 。 程 中 应 着 重 平衡 病床 资源 过 剩 以 及 病床 资源 紧张 的 问题 

分 配 问 题 。 安 佰 玲 等 中 根据 某 眼科 医院 的 病床 安排 指标 , 采 ”Wang 等 建立 以 医院 拒绝 病人 就 诊 次 数 最 少 和 等 待 时 所 花 成 
用 了 休克 疗法 的 思想 建立 非 线性 模型 ， 在 对 指标 进行 评分 后 ”本 最 少 为 目标 的 非 线 性 规划 模型 并 采用 元 模型 仿真 优化 [9。 
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徐 雷 等 在 考虑 病人 等 待 和 转 诊 成 本 的 前 提 下 ， 根 据 患 者 信息 
建立 混合 整数 线性 规划 的 病床 分 配 调度 模型 403。 宋 鸿 芳 等 考 
虑 患者 不 能 准时 入 院 或 提前 出 院 两 种 因素 ， 构 建 排队 指标 多 


周 晓 鸣 等 在 考虑 临时 床位 与 固定 床位 的 情况 下 提出 一 个 两 阶 
段 模型 , 采用 遗传 算法 求解 联合 分 配 问题 07。Bhattacharjee 和 
Raybg 提 出 对 病床 资源 分 配 和 患者 住院 时 间 相 结合 的 联合 优 
此 能 更 好 完善 患者 排队 系统 。 上 述 研究 主要 讨论 的 是 关于 病 
床 资源 的 多 目标 分 配 优化 方法 ， 能 较 好 地 刻画 患者 以 及 医院 
的 实际 情况 。 在 实际 情况 下 要 考虑 的 因素 更 为 复杂 ， 多 目标 
模型 的 关键 之 处 在 于 所 考虑 的 因素 能 和 否 准确 的 描述 患者 接受 
医疗 服务 的 过 程 。 目 前 研究 主要 考虑 患者 的 因素 较 多 ， 优 化 
的 目标 集中 在 降低 运营 成 本 以 及 患者 的 被 拒绝 就 诊 方面 ， 对 
于 院 方 与 患者 之 间 的 关系 ， 医 院 与 患者 之 间 的 公平 性 以 及 从 
医院 方 因素 考虑 的 研究 较 少 ， 有 和 较 大 的 研究 空间 。 
综 上 所 述 ， 本 文 着 重 在 考虑 病 患 公平 性 以 及 降低 成 本 的 
角度 进行 多 目标 的 分 配 优化 研究 。 在 体现 患者 公平 性 方面 ， 
本 文 根 据 相 关 的 医院 卫生 政策 9 纠 以 及 医院 运营 管理 政策 3 
创新 性 地 引入 一 个 权重 测度 指标 作为 公平 性 的 判断 标准 。 从 
院 方 和 角度， 本文 选取 医院 的 运作 成 本 作为 保障 医院 日 常 运作 
的 因素 。 在 模型 的 建立 过 程 中 ， 除 了 设置 运作 成 本 以 及 公平 
权重 最 优化 两 个 目标 函数 ， 结 合 病房 的 总 容量 ， 病 人 的 诊疗 
费用 以 及 系统 的 稳定 性 等 约束 条 件 构建 一 个 多 目标 随机 规划 
模型 MOP(multi-objective programming)， 并 通过 算 例 分 析 验 
证 了 模型 在 实际 运用 过 程 的 有 效 性 。 研 究 结 果 不 仅 可 以 拓展 
多 目标 下 的 医院 病床 分 配 研究 ,该 模型 的 合理 运用 有 助 于 医 
院 管理 者 在 实际 医院 管理 中 对 病床 资源 进行 合理 分 配 。 


1 问题 描述 


本 文 的 问题 考虑 病 患 的 公平 性 以 及 医院 的 运作 成 本 两 个 
天 素 ， 对 医院 单 科室 内 的 特殊 病人 与 普通 病人 的 病床 进行 分 
配 的 问题 。 病 人 根据 诊疗 要 求 的 诊疗 服务 分 为 普通 病人 和 特 
殊 病 人 。 医 院 住院 部 内 的 病房 依照 病人 需求 类 型 进一步 分 类 
为 公共 病房 与 私人 病房 。 考 虑 到 病房 分 配 是 一 个 较为 长 期 的 
过 程 ， 需 要 在 一 个 有 界定 的 时 间 范 围 内 解决 ， 因 此 限定 医院 
的 病房 类 型 在 半年 至 一 年 内 的 时 间 内 该 房间 类 型 不 会 改变 。 

医院 日 常 运作 成 本 主要 体现 在 医院 工作 人 员 的 开支 以 及 
医疗 设备 购买 维护 上 ， 为 了 方便 模型 的 建立 ， 本 文 汇 总 医院 
运作 产生 的 成 本 费用 ， 再 针对 每 种 不 同类 型 的 病人 进行 成 本 
全 
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平均 化 处 理 ， 即 服务 单位 病人 所 需 的 成 本 。 在 反映 反映 公平 
性 方面 ， 不 同 病 人 的 病情 具有 复杂 性 、 不 确定 性 ， 不 能 很 好 
地 度量 来 体现 医院 的 服务 质量 ， 因 此 本 文 不 考虑 是 否 要 预 留 
一 定 的 床位 来 保证 急诊 病人 的 病床 供给 ， 并 依据 此 设立 一 个 


部 门 间 的 病床 分 配 模型 ， 利 用 动态 规划 的 方法 取得 最 优 解 09。 


虑 病房 里 的 病床 容量 的 具体 份额 。 只 要 任何 类 型 病房 的 空间 
里 至 少 有 一 张 床位 ， 病 人 即 可 入 住 。 不 考虑 相同 性 别 的 患者 
一 定 要 住 在 同一 间 病 房 。 
表 1 模型 参数 说 明 

Tab.1 Model parameters description 


符号 说 明 
7 病人 类 型 集合 ， 用 指数 ;表示 
J 病房 类 型 集合 ， 用 指数 j 表示 
7 病人 完整 诊疗 下 来 的 所 需要 的 时 间 周 期 ， 用 指数 表示 


i=1,2 1 表示 VIP 病人 ，2 表示 普通 病人 
j=1,2 ”1 表示 私人 病房 ，2 表示 公共 病房 
i;iel,jeJ 表示 该 患者 为 i 型 患者 ， 确 诊 类 型 为 i 型 ， 住 在 j 型 病房 
tel[0,7] 完整 诊疗 的 时 间 周 期 
No 病房 的 总 数 
b, 了 类 型 病房 的 病床 数 
8j(*) 站 型 病 患 间隔 到 达 时 间 的 概率 分 布 
5 六 型 患者 的 日 均 到 达 率 
访 (z) 消 型 患者 住院 时 间 的 概率 分 布 
My 型 患者 平均 受到 服务 的 概率 
医院 服务 一 名 患者 的 单位 成 本 
人 立 型 患者 的 容忍 赣 值 
0 模型 中 的 随机 变量 
N; i 类 型 的 病人 分 配 到 的 7 型 病房 的 数量 
K; 型 病人 的 准 入 性 指标 
Ys j 型 病 患 的 预计 到 达 率 
Cy 型 病 患 的 响应 度 指 标 
D,(1,9) ”经 过 1 后 诊疗 完成 的 六 型 病 患 数量 
Wj(9) 六 型 患者 等 待 诊疗 的 时 间 
m=0,4.…,M 病人 预计 到 达 率 的 离散 化 指数 
M 病人 预计 到 达 率 的 离散 度 
hi 型 患者 的 离散 到 达 率 
1 表示 六 型 患者 的 计划 到 达 率 为 ho ， 
0 表示 计划 到 达 率 不 为 各 
1=01.…No 病房 数量 的 离散 化 指数 
Nm e{0,1} 1 表示 病房 的 数量 为 n，0 表示 病房 数量 不 为 7 
当 Kj 和 Nw 都 为 1 时 为 1， 
否则 为 0 
6 最 小 化 权重 指数 


3 ”模型 建立 
3.1 权重 测度 指标 

本 文 设立 了 结合 准 入 性 公平 以 及 响应 性 公平 的 权重 测度 
指标 来 衡量 病床 分 配 的 公平 性 。 根 据 医疗 政策 方面 0220 和 医 


权重 测度 指标 来 衡量 病床 分 配 公平 性 。 权 重 测度 指标 主要 从 
两 方面 来 考虑 : 病人 入 院 等 待 时 间 以 及 患者 所 能 接受 最 大 等 
待 天 数 。 
问题 的 主要 难度 在 于 如 何平 衡 医院 的 运作 成 本 与 患者 公 
平 性 ， 因 此 需要 建立 一 个 多 目标 随机 规划 模型 平衡 成 本 与 公 
平 ， 并 采取 数值 检验 改 模型 的 有 效 性 以 及 实用 性 。 
2 ”符号 说 明和 假设 
2.1 符号 说 明 

模型 参数 说 明 如 表 1 所 示 。 
2.2 ”模型 假设 

假设 Al: 不 同类 型 的 患者 相对 独立 , 同一 类 型 的 患者 所 
能 接受 最 大 等 待 天 数 相 同 ， 按 照 先 到 先 得 的 原则 接受 医院 诊 
疗 。 病 人 因 个 人 原因 的 延误 治疗 以 及 中 途 退 出 等 因素 不 予 考虑 。 

假设 A2: 每 种 类 型 的 病人 只 会 分 配 到 对 应 的 病房 。 不 考 


院 运营 策略 考量 2223， 设 立 一 个 指数 玉 作为 准 入 性 公平 的 度 
量 指标 ， 即 入 院 患者 确诊 为 住院 患者 的 比例 。 此 外 设 一 个 病 
人 所 能 接受 的 最 大 等 待 天 数 reR+ ， 不 确定 函数 变量 为 0 ， 等 
待 的 时 间 为 W( 90)，CVaR, (the conditional value-at-risk)24 为 条 
牛 风险 值 , 显示 了 医院 院 患者 确诊 为 住院 患者 的 转换 转换 率 ， 
其 次 根据 病人 的 等 待 时 间 分 布 以 及 所 接受 的 最 大 等 待 天 数 设 
立 指数 x 作为 响应 性 公平 度量 指标 ， 定 义 如 下 : 
定义 1 当 满 足 
CVaR, (0)=E[W(0)| W(0)>g,], 
K(W(0)<gq,)=a,0 el[0,1] 
Ssup{a <0| CVaR, <7} 
wa: 用 作 衡 量 病 人 等 待 诊疗 时 间 , 指 的 是 在 & 右 侧 检验 条 
件 分 布 中 的 预期 等 待 时 间 , 指标 a.(9) 是 指 满足 x 的 右 侧 检验 
下 病人 预计 等 待 诊疗 时 间 不 超过 病人 的 容忍 癜 值 的 最 大 值 。 
该 指标 从 预计 等 待 时 间 这 个 角度 反映 了 医疗 服务 的 水 平 。 指 
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标 a.(9) 有 如 下 的 性 质 。 
a) 假 设 W(0)>W(80)， 则 a.(80)<a(0); 
b) 假 设 W(9)<r+ ， 则 @.(9)=1; 
co) 假设 BW(9)]>r ， 则 (9)=0; 
d) 0.(0)< K(W(0)<7)。 
证 明 1) 假设 W(0)>W(0)， 
CVaR, (0)>CVaR,(0,), 则 a(0)s0(0); 
2) 假 设 W(O)srz ， 则 有 CVaR(0 )<r ， 那 么 w (9O)=1; 
3) 假 设 EIw(o)jlj>r ， 则 有 CyaR(g)=E[W(9)]>r ， 那 么 


有 任意 ae[0,1] ， 
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接收 率 。 式 (8) 表 示 医 院 接 收 每 一 类 病人 的 诊疗 程度 由 患者 的 
日 均 到 达 率 和 预计 到 达 率 的 除 积 体现 。 


Ay 
<lviel,jeJ 
bj;Ni ns 7 (7) 


院 运作 成 本 的 病床 配置 优化 研究 


Vi ” 
Ky = viel,jeJ (8) 


也 


式 (9) 为 式 (1) 到 式 (8) 的 所 有 变量 的 取 值 范围 。 

Ni e N;K;,,0; e R;0 < Ki,o Slviel,jeJ] (9) 

模型 中 虽然 考虑 了 病房 的 数量 ， 但 是 病房 有 时 候 需 要 进 

行 扩容 或 者 翻新 等 处 理 ， 式 (3) 也 可 以 由 每 个 房间 所 需要 的 日 


cc(9)=0; 

4) 令 由 =K(W(g)<rz)，cyam(g)=EIY(o)Iw(9)>zj>z ， 则 有 
a (0)<K(W(O)<1) 。 

定义 2 结合 准 入 性 指标 以 及 响应 性 的 指标 & ， 将 公 
平 性 的 指标 定义 为 K, ，K, =Kxa 。 


准 入 性 指标 大 和 响应 性 指标 2% 是 分 别 用 来 衡 0 


诊 服务 的 两 个 连续 过 程 ， 两 者 相 乘 可 得 一 个 衡量 公平 性 的 权 
重 测度 指标 。K。 权重 测度 指标 表现 og ed 
预期 等 待 确诊 时 间 的 关系 ， 为 保障 多 个 目标 之 间 的 平衡 ， 将 
K。 作为 公平 最 大 化 衡量 标准 。 


3.2 多 目标 随机 规划 模型 
3.2.1 目标 函数 

式 (]) 表 示 了 医院 的 运作 成 本 最 低 的 目标 函数 。 目 标 函 数 
表示 了 在 预期 时 间 窗 T 内 完成 诊疗 的 所 有 类 型 患者 之 和 
Dj(7,0) 与 诊疗 单位 病人 的 成 本 费用 广 相 乘 的 总 和 。 


min f =22 ElD, (7,0)] (1) 


于 本 文 要 达到 公平 性 最 优 ， 所 以 取 最 大 值 来 表示 权重 
测度 指标 可 提高 的 最 大 限度 ， 以 最 大 程度 的 提升 医院 的 公平 
性 。&w 如 3.1 小 节 所 述 ， 分 别 代 表 了 权重 测度 指标 Ko 
的 准 入 性 公平 指数 与 响应 性 公平 指数 ， 即 病人 的 就 诊 率 以 
及 集合 病人 等 待 时 间 分 布 的 病人 等 待 闷 值 。 式 (2) 描 述 了 公平 
E 权 重 测度 指标 最 优 的 目标 ，Ks% 为 准 入 性 公平 指数 与 响应 
性 公平 指数 的 乘积 


重 疯 


= 


oa O) 
3.2.2 约束 条 件 
(1) 式 (3) 表 示 了 关于 病房 总 容量 的 约束 ， 为 了 确保 分 配 
的 病房 的 容量 NV 不 超过 病房 的 总 数 No 。 


Zo <No G) 


(C) 式 (4) 表 达 的 是 诊疗 时 间 工 结束 前 完成 治疗 的 立 型 患 
者 的 数量 。 由 于 患者 的 计划 到 达 率 与 住院 的 时 间 分 布 没有 规 
律 性 ， 所 以 很 难 采用 有 具体 的 封闭 函数 表达 式 来 描述 这 种 函数 
关系 ， 暂 用 函数 2 来 表示 这 段 闭 合 函数 关系 式 。 在 接 下 来 的 
线性 化 求解 中 进行 进一步 的 求解 。 

D;,(T,0)=Q(N;,,7s, 8; (x),h(x)) Viel, je (4) 

(G3) 式 (5) 和 式 (6) 是 3.1 小 节 所 提 到 过 的 公平 性 权重 测度 
指标 。 用 简写 C 来 代替 CVaR, 。 设 立 病人 所 能 接受 的 最 大 等 
待 天 数 teRi 为 一 个 容忍 闵 值 , 不 确定 函数 变量 为 9, 等 待 的 
时 间 为 W(0)。 式 (5) 表 示 了 通过 检验 后 的 响应 性 指标 。 式 (6) 
中 右 侧 的 变量 内 容 分 别 为 六 型 病 患 间隔 到 达 时 间 , 住院 时 间 ， 
预计 到 达 率 分 布 ， 以 及 i 类 型 的 病人 分 配 到 的 j 型 病房 的 数 
量 ， 表 示 了 患者 的 确诊 率 。 

0 =sup{a >0| C(W, (0))<r,}viel, je (5) 
C(W,(0))=C(N;,7y,8s (x),h, (x))viel, je (6) 
3.2.3 系统 的 稳定 性 约束 

为 了 保持 整个 多 目标 随机 规划 模型 以 及 整个 医院 诊疗 流 
程 的 稳定 性 ， 需 要 保证 式 (7) 的 服务 强度 要 小 于 1.0， 即 六 型 
病 患 的 服务 接收 率 要 远 远 小 于 在 公共 病房 中 所 有 患者 类 型 的 


常 运营 成 本 来 呈现 ， 但 是 在 本 文中 不 予 考虑 。 
3.3 线性 化 转换 

式 (4) 和 式 (6) 缺 乏 确定 的 函数 表达 式 , 所 以 本 文 根 据 模型 
之 间 各 个 变量 的 关系 使 用 线性 化 方法 将 多 目标 随机 规划 模型 
转变 成 混合 整数 线性 规划 模型 。 当 六 型 病 患 的 预计 到 达 率 方 
和 i 型 病人 住 在 j 型 病房 的 数量 入; 车 人 出， Zi(7,9) 和 
CVaR,(W, (9)) 的 关系 ， 如 式 (10) 所 示 。 进 一 步 对 六 型 病 患 的 预 
计 到 达 率 进行 分 解 ， OU 0-1 变量 表示 Ns 。 


= -Sr hn An = et; 
jee ,M} (10) 
Ni =>, 2 = 上 
Ni e{0,1l} viel, jeJ,ne{0,,...,No} (11) 


在 上 式 中 参数 各 表示 的 是 立 型 患者 的 离散 到 达 率 , 而 参 
数 Km 和 Nw 的 组 合 都 代表 着 病房 病床 的 容量 和 病人 到 达 率 的 
组 合 。 通 过 前 文 参数 来 对 应 的 计算 [Dow(7,9)] 和 wm。 令 
D5,=B[Dym(7.9)| ,将 上 述 的 参数 关系 之 间 的 线性 化 转换 如 式 
(12)- 式 (21) 所 示 。 只 要 病人 预计 到 达 率 的 离散 程度 M 够 大 ， 
则 该 模型 的 解 近似 于 原来 多 目标 随机 规划 模型 的 解 。 


M No 


min fi= > > > > pi Xi Dom (12) 
iel jeJm=0n=0 
maxfy =& (13) 
No 
D2 Nons No (14) 
iel je/ n=0 
6<$ itm 六 ar nVviel,jeJ (15) 
力 > 各 Ky, < Mb; SN vi el,jeJ (16) 
m=0 n=0 
Xmn < KmVieL, jeJ,net{0,l,...,No},me{0,l,...,M} (17) 
No, Mo, 
Kyn =1, 2 Ny =1, > xinn =1Viel, jsy (18) 
m=0 n=0 m=0n=0 


Nm > Kom + No ~1Viel, je Tne{0 lb Nome{0,l,...,M} (19) 
G eR,0<E <1;Ki, Nin, Xomm {0,1}Viel, jeJ,net{0,l,...,No} (20) 
me{0,1,...,M} (21) 

4 ”算法 设计 
为 求解 3.3 小 节 的 多 目标 整数 线性 化 模型 ， 下 1 多 目标 
算法 是 一 种 较 好 的 选择 。 其 中 , NSGA2 算法 被 广泛 运用 于 复 
杂 的 ， 非 线性 的 多 目标 优化 问题 5221。 但 NSGA2 在 遗传 各 
群 迭 代 的 收敛 性 ， 运 算 效 率 ， 防 止 陷 入 局 部 最 优等 方面 还 需 
改进 。 本 文通 过 在 快速 非 支 配 排序 中 结合 拥挤 度 信 息 和 种 群 
内 部 分 布 性 保持 方法 两 个 因素 定义 新 的 个 体 快速 支配 强度 排 
序 赋值 ， 在 精英 保留 策略 的 引入 伪 适 应 度 ， 在 交叉 变异 过 程 
中 根据 种 群 个 体 的 排序 赋值 信息 重新 改进 交叉 算 子 等 方面 对 
原 算法 进行 改进 。 
4.1 编码 
根据 医院 单 科室 病床 分 配 的 特点 ， 本 文 在 编码 中 选择 病 
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组 ，X=[A,B]，A={aeAl<si<N,0<a<1}B={h eBl<i<N,0<b <1) 。 
4 表示 的 是 患者 的 到 达 率 , 即 一 组 1 XN 的 矩阵 , NN 为 到 
达 医 院 的 患者 数量 ， a 为 表示 i 个 患者 的 分 配 优 先 级 。B 表 
示 每 个 患者 在 接受 诊疗 之 后 分 配 到 各 个 病房 的 集合 ， 也 是 一 
组 1XN 大 小 的 矩阵 ，N 为 接受 诊疗 后 的 患者 数量 ， 户 为 判 
断 第 i 个 患者 分 配 到 私人 病房 或 公共 病房 。 
4.2 改进 的 快速 支配 排序 
在 NSGA2 中 采用 的 是 快速 非 支 配 排序 法 ， 但 在 求解 高 
维 复杂 的 混合 整数 规划 模型 时 由 于 多 目标 的 数量 众多 导致 排 
序 时 间 过 长 ， 且 非 支 配 解 集 中 存在 大 量 的 无 效 非 支 配 解 降低 
了 种 群 分 层 排序 的 收敛 效率 。 本 文采 用 改进 的 快速 支配 排序 ， 
将 原 有 的 比较 个 体 数 增加 为 两 个 ， 并 引入 新 的 支配 判别 强度 
指标 9 判别 种 群 的 优越 性 。 指 标 计算 公式 如 下 : 


6- 让 2 和 C2) 


2 
其 中 N 为 目标 的 个 数 ， je 和 fr 分 别 表示 解 集中 第 个 
标的 最 大 与 最 小 值 , 通过 式 (21) 排 除 解 集 中 的 无 效 支配 解 , 在 
非 支配 解 集中 单个 个 体 的 支配 强度 越 小 ， 则 优化 效果 越 好 ， 
确保 优异 个 体能 进入 遗传 迭代 。 
4.3 改进 的 拥挤 度 计 算 
NSGA2 算法 采用 计算 拥挤 度 来 保障 种 群 的 多 样 性 , 一 般 
将 种 群 中 某 个 点 的 周围 定点 的 密度 称 作 拥挤 度 281。 虽 然 拥挤 
度 能 保持 种 群 的 均匀 分 布 ， 但 如 图 1 所 示 ， 个 体 a 和 个 体 b 
之 间 的 距离 相 较 于 a 和 另 一 点 的 距离 更 短 ，a 和 的 拥挤 
度 距离 计算 较 大 有 旦 相似 ， 可 能 会 同时 保留 或 淘汰 。 该 策略 会 
导致 分 布 较 好 的 某 些 个 体 淘 汰 ， 而 分 布 较 差 的 个 体 保留 。 


图 1 拥挤 度 距 离 比较 


Fig. 1 Crowding distance comparison 


为 了 保持 种 族 的 分 布 性 ， 采 | 
挤 距离 计算 方法 ， 公 式 如 下 : 


He 2 | 


Df — pa 
lf- 六 A 于 
n 


中 分 子 到 j- 用 | 为 原 有 拥挤 度 计算 公式 ， /和 大 


分 别 表示 个 体 i+1; 六 7 在 对 目标 排序 后 在 第 个 子 目标 上 的 
目标 函数 值 。 

.4 改进 的 精英 保留 策略 

在 NSGA2 算法 中 ， 精 英 保留 策略 依照 非 支 配 关 系 进行 
排序 ， 支 配 等 级 低 的 个 体 填充 到 新 的 父 代 集 合 中 ， 直 到 种 群 
数量 达到 N， 较 大 程度 的 保留 进化 过 程 优秀 个 体 并 防止 其 流 
失 。 但 由 于 传统 的 精英 策略 保留 是 固定 值 ， 非 支配 解 集 收敛 
效果 不 佳 ,为 了 提高 解 集 的 收敛 效率 , 引入 伪 适 应 度 值 09， 


种 考虑 拥挤 度 方差 的 拥 


5= 


(23) 


| 中 


A 


上 


全 
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人 的 计划 到 达 率 以 及 医院 科室 的 病床 分 配 能 力作 为 两 个 决策 。 增加 精英 的 保留 规模 ,f 越 小 则 说 明 个 体 越 优秀 , 计算 公式 如 下 : 
变量 进行 编码 。 病 床 的 容量 以 及 病人 的 离散 化 计划 到 达 率 既 vo; 

可 用 整数 表示 ， 也 可 用 编码 的 方式 呈现 。 因 为 NSGA2 算法 Fok+tl_ 和 于 (24) 
无 法 直接 对 问题 空间 内 的 参数 进行 处 理 ， 所 以 本 文 对 连续 性 

变量 进行 离散 化 后 的 离散 型 变量 浮 点 数 编码 。 染 色 体 分 成 两 其 中 :为 帕 累 托 前 沿 的 阶 数 ;5 为 拥挤 度 距 离 ; N 为 目标 函 


数 的 数量 。 通 过 改进 ， 在 种 群 迭代 前 期 非 支 配 解 集 较 少 时 ， 
较 好 的 保留 精英 个 体 ， 在 后 期 非 支 配 解 集 较 多 的 情况 下 能 扩 
大 精英 个 体 保留 的 规模 ， 使 得 种 群 快速 收敛 。 

4.5 交叉 和 变异 
由 于 传统 的 SBX(Simulated Binary Crossover) 交 叉 算 子 的 
全 局 搜索 能 力 较 弱 BU， 为 了 更 好 地 保证 种 群 的 多 样 性 ， 采 用 
算术 交叉 分 子 蔡 代 传统 的 SBX 交叉 算 子 ， 计 算 公 式 如 下 ， 


Brank 


人 Arank + Brank 25) 
% =05[d4+o)s +(—a)5,] (26) 
2 =0.5[(—a)S, +(+oa)S,] (27) 


其 中 :Arank 和 Brank 分 别 表示 个 体 4 与 个 体 B 的 非 支 配 排 
序 值 ; 交叉 算 子 系数 a 表示 每 个 个 体 在 快速 支配 排序 强度 中 
的 排序 关联 度 。 在 运算 前 期 & 的 变化 较 大 ， 非 支配 排序 值 较 
小 的 个 体 占据 较 大 数量 ， 随 着 运算 的 推进 ， 群 内 的 个 体 逐 渐 
收敛 于 同一 帕 累 托 前 沿 ， cx 趋 近 于 0.5。 式 (25) 和 式 (26) 选 取 


父 代 个 体 5,、5, 交叉 得 到 子 代 5 、 > 。 
变异 操作 则 从 父 代 种 群 中 随机 选取 父 代 个 体 5,， 产 生子 
人 ws ， 多 项 式 变异 公式 如 下 ， 
$3=53+[Q, 一 1]( xu — Xin ) (28) 
其 中 : % 为 变异 分 布 指数 ; xs 和 xm 为 变量 的 上 下 限 。 


4.6 算法 步骤 
步骤 1: 对 所 有 参数 初始 化 ， 随 机 生成 含有 V 个 个 体 的 
种 群 S$， 任 意 选 定 一 个 个 体 作 最 优 解 ， 记 为 b。 
步骤 2: 引入 快速 支配 强度 排序 法 对 父 代 种 群 5 的 内 部 
个 体 排序 ， 根 据 排序 支配 等 级 选取 个 体 保留 到 外 部 集合 。 
步骤 3: 对 种 群 5, 中 个 体 进行 改进 的 拥挤 度 计算 ， 赋 予 
个 体 相应 的 非 支 配 等 级 计算 个 体 拥挤 度 。 
步骤 4: 运用 伪 适 应 度 值 的 精英 保留 策略 对 种 群 5, 进行 
排序 ， 对 种 群 5, 交叉 变异 后 得 到 包含 N 个 个 体 的 新 子 种 群 ， 


记 为 0， 将 这 两 个 种 群 合 并 ， 得 到 个 体 数 为 原来 2 倍 的 种 群 
5S ， 挑 选 出 前 个 个 体 组 成 第 S+17 代 种 群 ， 并 将 其 产生 的 最 
优 个 体 记 为 b1。 

步骤 5$: 将 bl 与 b 进行 比较 , 车 bl 优 于 个 体 b， 则 选用 
个 体 bl 。 


步骤 6: 判断 是 否 达到 最 大 友 代 次 数 ， 如 果 是 ， 
步骤 7; 否则 执行 步骤 2。 

步骤 7: 输出 个 体 最 优 解 。 

算法 步 又 流程 图 如 图 2 所 示 。 
5 ”算法 测试 和 算 例 分 析 
5.1 多 目标 算法 性 能 评价 指标 

为 验证 改进 后 的 算法 性 能 , 本 文选 用 Zitzler 和 Deb 提出 
BJ 中 的 ZDT1、ZDT2、ZDT3 和 ZDT6 测试 函数 进行 实验 ， 
并 与 NSGA2 测试 结果 进行 比较 。 同 时 为 了 检验 改进 后 的 算 
法 在 求解 混合 整数 线性 问题 时 的 全 局 搜索 优化 和 约束 处 理 能 
力 , 选用 DTLZ1、DTLZ4、DTLZ6、DTLZ7 函数 测试 NSGA2- 
DS 在 求解 过 程 的 收敛 性 与 多 样 性 B43。 在 评价 NSGA2-DS 方 
选取 三 个 指标 来 分 析 改 进 后 算法 的 收敛 性 、 多 样 性 、 以 及 
综合 性 ， 指 标 如 下 所 示 ， 
(1) 收 敛 性 的 评价 指标 世代 距离 (GD) 
GD 表示 解 集 到 真实 的 帕 累 托 前 沿 上 最 近 参 考点 的 最 小 
欧式 距离 中 以 评估 算法 的 收敛 性 ， 计 算 公 式 如 下 ， 


则 执行 
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8 =mipF (x)-F(p),x es (29) 
癌 
Bi 
GDGS.P)=- HT G30) 


其 中 : 8 表示 的 是 x 与 帕 累 托 近 似 前 沿 P 了 上 的 距离 最 近 的 参 
考点 p 之 间 的 欧式 距离 ;5 为 算法 所 求 得 的 解 集 ， 世 代 距 离 
值 越 小 表示 算法 收敛 性 越 好 。 
CO) 多 样 性 的 评价 指标 空间 指标 (SP) 
选用 Schott JR64 提 出 的 空间 指标 度量 解 集 内 各 个 解 到 
其 他 解 之 间 最 小 距离 的 标准 差 来 衡量 解 在 空间 的 分 布 情况 ， 
计算 公式 如 下 ， 


gi -mil 多 (x')—F Co))* ES,xizx (31) 


，. 冉 s 2 
Spacing(P)=, i —g;) (32) 


8 表示 的 是 x' 到 ww 这 两 个 连续 解 之 间 的 曼哈顿 距离 ，&g 
表示 所 有 8 的 均值 ，SP 为 第 8; 个 解 到 其 他 的 解 的 最 小 距离 ， 
SP 的 值 越 小 ， 解 集 的 多 样 性 越 好 。 
(G3) 综 合 性 的 评价 指标 改进 后 的 反 向 迭代 距离 LIGD7) 

反 向 迭代 距离 B51(Inverted Generational Distance) 表 示 帕 
累 托 前 沿 上 的 每 个 参考 点 到 解 集中 邻 域内 最 近 的 解 的 均值 距 
离 ， 在 计算 过 程 中 当 参 考点 较 少 时 ， 如 图 3 所 示 ， 解 集 9 中 
的 所 有 解 都 受到 解 集 % 的 支配 并 且 解 集 % 的 质量 明显 高 于 
解 集 $s， 如果 采用 IGD 评价 会 得 到 解 集 5; 的 质量 高 于 解 集 
51 的 情况 。 


图 2 改进 后 的 NSGA2-DS 算法 流程 图 
Fig.2 Flow chart ofthe improved NSGA2-DS algorithm 
因此 本 文 引入 改进 后 的 IGD*B9I， 计 算 公 式 如 下 ， 


G33) 


d+ 
1GD:(P,S)=—— 34 
IP| ( ) 


IGD' 在 原 有 IGD 计算 参考 点 p 到 解 x 的 距离 的 同时 考 
虚 支 配 关 系 ，d: 是 考虑 解 的 支配 关系 的 帕 累 托 前 沿 上 参考 点 
p 到 与 之 最 近 的 解 的 欧式 距离 ，IGD' 在 缺少 帕 累 托 前 沿 参考 
集 的 情况 下 也 一 样 适用 ， 在 实际 优化 过 程 中 的 适用 性 更 高 。 
IGD'* 的 值 越 小 ， 说 明 算 法 的 综合 性 能 更 好 。 
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上 
图 3 IGD 的 局 限 性 
Fig.3 Limitations of IGD 


5.2 算法 测试 
5.2.1 测试 参数 设置 

为 检验 改进 后 的 算法 在 收敛 性 、 多 样 性 、 分 布 性 等 方面 
的 性 能 ， 将 算法 与 原来 的 NSGA2 进行 对 比 。 对 两 种 算法 均 
做 如 下 实验 参数 设置 : 种 群 规模 N 设置 为 100， 种 群 交叉 分 
布 指数 a.=1， 交 叉 参 数 设置 为 20， 交 又 概 率 为 0.9， 种 群 变 
异 分 布 指 数 mw =1， 变 异 参数 设置 为 20， 变 异 概率 为 0.1， 最 
大 迭代 次 数 maxFE 为 500 次 。 在 NSGA2-DS 算法 中 伪 适 应 
度 值 设置 为 fs[0.2.0.8] ， 保 证 精英 个 体 的 增长 规模 。 对 ZDT 
系列 测试 函数 中 的 ZDT1 和 ZDT2 的 变量 个 数 设 置 为 30， 
ZDT3 与 ZDT6 的 变量 个 数 设 置 为 10，DTLZ 系列 测试 函数 
的 变量 目标 设置 为 5， 每 组 实验 测试 独立 运行 30 次 。 
5.2.2 实验 测试 结果 及 分 析 

两 种 算法 在 ZDT 系列 测试 函数 的 数值 仿真 结果 如 图 4~7 
所 示 ，DTLZ 系列 测试 函数 的 数值 仿真 结果 如 图 8~11 所 示 。 
由 于 ZDT 系列 函数 的 决策 变量 数目 较 多 , 适合 比较 高 维 的 决 
策 变量 的 混合 整数 线性 规划 问题 。 由 图 4~7 可 知 ，ZDT1、 
ZDT2、ZDT3 测试 函数 的 求解 结果 显示 NSGA2-DS 算法 和 
NSGA2 算法 都 能 得 到 较 好 的 收敛 效果 , 但 是 改进 的 算法 明显 
比 NSGA2 收敛 结果 更 接近 帕 累 托 前 沿 ， 在 求解 ZDT6 测试 
函数 时 ，NSGA2-DS 的 收敛 效果 明显 优 于 NSGA2 算法 且 
NSGA2 求解 结果 中 存在 较 多 非 支 配 解 。 


隔 


0 02 0.4 0.6 08 1 12 
图 4 NSGA2-DS 与 NSGA2 的 ZDTI1 测试 函数 对 比 
Fig.4 Comparison ofthe ZDT1 test function of NSGA2-DS and NSGA2 
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图 5 NSGA2-DS 与 NSGA2 的 ZDT2 测试 函数 对 比 
Fig.5 Comparison ofthe ZDT2 test function of NSGA2-DS and NSGA2 
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图 6 NSGA2-DS 与 NSGA2 的 ZDT3 测试 函数 对 比 
Fig.6 Comparison ofthe ZDT3 test function of NSGA2-DS and NSGA2 


图 11 NSGA2-DS 与 NSGA2 的 DTLZ7 测试 函数 对 比 
|。 Fig. ll Comparison ofDTLZ7 test functions of NSGA2-DS and NSGA2 


函数 的 GD、SP、IGD-* 指 标 进 行 统计 ， 每 个 测试 函数 都 运行 
oa * 30 组 的 运算 。 如 表 5 所 示 ， 本 文 给 出 了 部 分 NSGA2-DS 算 
02| 法 在 计算 IGD' 指 标的 实验 数据 ， 表 2- 表 4 分 别 统 计 两 个 测 
"| 试 函数 的 GD、SP、IGD+ 指 标 均 值 与 均 方 差 。 由 表 2 可 知 ， 
EY | e 除了 ZDT3，ZDT6，DTLZ7 的 NSGA2-DS 均值 大 于 NSGA2 
图 7 NSGA2-DS 与 NSGA2 的 ZDT6 测试 函数 对 比 的 求解 结果 , NSGA2-DS 的 收敛 性 结果 总 体 由 于 NSGA2; 均 
Fig. 7 Comparison of the ZDT6 test function of NSGA2-DS and NSGA2 方差 中 除了 DTLZ4 和 DTLZ6，NSGA2-DS 的 结果 明显 小 于 
如 图 8~11 所 示 ， 渐 变色 块 为 DTLZ 系列 测试 函数 的 真 NSGA2， 有 具有 更 强 的 鲁 棒 性 。 同 理 ， 由 表 3 可 得 NSGA2-DS 
实 帕 累 托 边界 , NSGA2-DS 和 NSGA2 算法 在 求解 DTLZ4 六 解 的 多 样 性 与 分 布 性 总 体 优 于 NSGA2， 表 4 的 IGD+ 指 标 对 
试 函数 时 的 收敛 效果 较 好 ; 在 求解 DTLZ1、DTLZ7 测试 函数 比 中 可 看 出 NSGA2-DS 的 综合 收敛 性 、 多 样 性 、 延 展 性 均 优 
的 结果 中 ，NSGA2-DS 的 收敛 性 、 所 求解 分 布 均匀 性 、 多 样 于 NSGA2。 虽 然 改 进 后 的 算法 在 各 个 指标 上 的 表现 总 体 优 于 
性 均 优 于 NSGA2， 但 在 DTLZ6 测试 函数 的 所 求解 收敛 性 与 改进 前 , 但 是 对 比分 析 后 得 出 NSGA2-DS 由 于 引入 伪 适 应 度 
分 布 性 均 差 与 NSGA2 算法 ， 存 在 较 多 伪 支 配 解 不 能 收敛 。 值 的 改进 精英 选择 策略 在 求解 前 期 的 非 支 配 解 较 多 ， 精 英 个 
体 保 留 数量 较 大 导致 收敛 过 快 易 造成 局 部 最 优 ， 但 总 体 上 
NSGA2-DS 的 性 能 优 于 NSGA2， 在 求解 复杂 高 维 问题 上 明 
显 优 于 NSGA2。 
表 2 算法 的 GD 指标 对 比 
Tab.2 Comparison of GD metrics of algorithms 
均值 均 方 差 
NSGA2-DS NSGA2 NSGA2-DS NSGA2 
ZDT1 1.2905E-03 1.1155E-01 1.0383E-01 2.2383E-01 
ZDT2 7.7206E-01 8.0814E-01 2.1526E-01 4.7526E-01 
图 8 NSGA2-DS 与 NSGA2 的 DTLZ1 测试 函数 对 比 ZDT3 3.4035E-01 1.1853E-01 2.9865E-02 1.4833E-01 
Fig. 8 Comparison of DTLZ1 test function of NSGA2-DS and NSGA2 ZDT6 1.2527E+00 7.6714E-01 6.7582E-02 8.7535E-01 
DTLZ1 1.0644E-01 1.1650E-01 1.1753E-01 1.1960E-01 
DTLZ4 5.6410E-04 1.6629E-03 3.7653E-03 7.2703E-04 
DTLZ6 1.2322E-06 3.3224E-02 5.2175E-01 4.2736E-02 
DTLZ7 1.7736E-02 6.3030E-03 8.7163E-04 9.7526E-04 
表 3 算法 的 SP 指标 对 比 
Tab.3 Comparison of SP metrics of algorithms 
均值 均 方 差 
NSGA2-DS NSGA2 NSGA2-DS NSGA2 
ZDT1 1.1857E-02 1.4286E-02 4.8513E-03 5.1675E-03 
图 9 NSGA2-DS 与 NSGA2 的 DTLZ4 测试 函数 对 比 ZDT2 8.8445E-03 4.3118E-02 2.6579E-02 3.1536E-02 
Fig.9 Comparison of DTLZ4 test function of NSGA2-DS and NSGA2 ZDT3 9.6351E-03 9.6484E-02 6.2306E-02 2.0089E-01 
本 ZDT6 6.4859E-02 8.1663E-02 4.7976E-02 2.9267E-02 
DTLZ1 6.2130E-02 1.2700E-01 3.5439E-02 6.7673E-02 
DTLZ4 5.2798E-02 4.9750E-02 1.0575E-03 2.4318E-02 


| ew 为 了 定量 比较 算法 性 能 ， 均 对 两 个 算法 在 求解 系列 测试 
ww 


tt 
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引 


MOP 


MOP 


下 民 2 DTLZ6 1.0868E-02 2.8241E-02 2.1739E-02 3.2587E-02 
5 DS 的 by DTLZ7 7.8114E-02 6.5458E-02 7.3S16E-03 8.6576E-03 
- \ 一 5.3 算 例 结果 及 分 析 
一 一 pa A 态 本 文选 取 某 家 公立 医院 数据 作为 数据 来 源 ， 数 据 时 效 性 
| 为 一 个 月 。 采 用 NSGA2-DS 和 NSGA2 算法 求解 出 多 目标 整 
图 10 NSGA2-DS 与 NSGA2 的 DTLZ6 测试 函数 对 比 数 规划 模型 的 结果 。 另 外 ， 本 文 还 引入 吴涛 97 等 提出 的 改进 


Fig. 10 Comparison of DTLZ6 test function of NSGA2-DS and NSGA2 的 黑洞 进化 算法 (Multi Objective Black Hole) 对 算 例 进 行 求解 。 
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黑洞 算法 多 目标 进化 算法 是 一 种 基于 信息 和 自 适应 化 策略 的 表 4 算法 的 IGD+ 指 标 对 比 表 
粒子 群 算法 ， 在 收敛 速度 ， 种 群 的 多 样 性 ， 种 群 的 收敛 性 县 Tab.4 Comparison table of IGD+ metrics of algorithms 
有 较 好 的 性 能 ， 因 此 将 MOBH 算法 与 改进 后 的 NSGA2-DS MOR 均值 均 方差 
算法 的 求解 结果 对 比 验 证 算法 的 前 沿 性 以 及 有 效 性 。 考 虑 到 NSGA2-DS ”NSGA2 NSGA2-DS NSGA2 
匠 累 托 最 优 解 的 情况 ， 本 文中 设置 种 群 规 模 为 100， 交 叉 概 ZDT1 2.6380E-01 2.6731E-01 2.1354E-01 2.7446E-01 
率 为 0.8， 变 异 概率 为 0.2， 最 大 友 代 次 数 为 500。 本 文 所 做 ZDT2 7.2289E-02 1.6480E-01 2.3286E-01 1.9115E-01 
的 数值 实验 均 在 Inter(R)Core(TM)i5-7300HQ 处 理 器 上 进行 ， ZDT3 7.7883E-03 7.4604E-02 1.7964E-04 2.1827E-01 
CPU 为 2.50GHz。 ZDT6 4.0492E-03 2.4583E-01 7.6553E-02 3.5836E-01 

不 同类 型 的 病人 有 不 同类 型 的 住院 时 长 、 在 所 能 接受 DTLZ1 8.9427E-01 2.9019E-01 5.2123E-02 3.5009E-01 
的 最 大 等 待 天 数 、 到 达 率 ， 同 时 医院 服务 不 同 的 病人 所 需 DTLZ4 7.2021E-02 1.8170E-01 5.4791E-03 1.9456E-01 
的 成 本 也 不 同 。 部 分 病人 的 基本 信息 如 下 表 6 所 示 ， 病 房 DTLZ6 5.1624E-03 4.2612E-01 5.2175E-01 5.1147E-01 
总 数 为 300 个 。 DTLZ7 8.3260E-02 1.4033E-01 6.8030E-04 1.2387E-01 


表 5 部 分 NSGA2-DS 算法 的 IGD+ 指 标 实验 数据 
Tab.5 Experimental data of IGD+ metrics for some NSGA2-DS algorithms 


ZDTI1 ZDT2 ZDT3 ZDT6 DTLZ1 DTLZ4 DTLZ6 DTLZ7 
1 1.2369E+00 5.5009E-01 8.3731E-04 1.1201E+00 1.3793E+00 7.4869E-02 4.5272E-02 2.6386E-01 
2 1.4817E-01 1.8469E-01 1.3343E-03 7.7460E-01 4.0639E-01 5.4145E-01 6.1207E-02 2.1322E-01 
3 1.3494E-01 1.2456E-01 5.5377E-03 5.4134E-01 5.4101E-01 7.2493E-02 4.9200E-02 2.5220E-01 
4 9.3185E-02 2.0480E-01 5.1773E-02 1.0118E+00 4.7429E-01 6.6197E-02 1.2824E-02 2.4571E-01 
5 4.6381E-01 6.6873E-02 1.9210E-02 5.7931E-01 2.4599E-01 6.9159E-02 1.6820E-02 3.2278E-01 
6 5.6183E-01 5.4086E-02 1.5812E-02 1.9935E-03 5.9775E-02 1.1230E-01 1.7746E-03 8.2786E-02 
7 4.6979E-01 2.8209E-01 1.0978E-03 1.4068E-01 7.7372E-02 6.8710E-02 1.3956E+00 5.7351E-03 
8 4.8367E-01 6.8352E-02 5.2100E-03 9.2168E-02 6.9379E-02 1.3291E-01 1.4465E+00 3.4214E-02 
9 8.0304E-02 9.5271E-02 3.1011E-03 3.0865E-01 6.5800E-02 6.5710E-04 6.8387E-01 1.1199E-02 
10 1.1164E-01 1.6310E-01 1.2379E-03 3.9021E-02 5.4966E-02 1.7579E-03 1.1603E+00 7.2705E-03 
表 6 部 分 病人 类 型 与 所 住 病房 的 基本 信息 表 100 时 ， 问 题 的 解 并 发 生 没 有 太 多 的 变化 ， 因 此 在 NSGA2- 
Tab.6 Table ofbasic information of some patient types and wards DS 算法 中 将 M 设置 在 100， 图 12 也 显示 了 使 用 NSGA2 与 
病人 类 型 所 住 病房 每 日 到 达 率 住院 时 长 /天 诊疗 成 本 /元 ”容忍 闽 值 /天 MOBH 算法 下 的 所 求 得 的 帆 累 托 最 优 解 ， 求 解 结果 与 
普通 公共 4.383 9.325 231.52 14 NSGA2-DS 相 比 较 。 结 果 显 示 在 总 体 求解 结果 上 NSGA2-DS 
特殊 私人 0.442 7.662 724.35 7 算法 与 MOBH 算法 求解 结果 均 优 于 NSGA2。NSGA2-DS 算 
普通 。 公共 3.523 2.617 502.31 6 法 的 所 求 得 最 优 解 效果 最 好 , MOBH 算法 虽然 在 求解 所 耗 时 
普通 。” 私人 0.235 2.701 219.23 3 比 NSGA2-DS 短 且 较 快 收敛 ,但 求解 结果 比 NSGA2-DS 略 
普通 。 公共 ”8.965 9.970 193.25 6 差 ， 在 公平 权重 0.5 相 邻 段 时 运作 成 本 相 较 于 NSGA2-DS 算 
村 殊 公共 0.178 10259 595.00 4 法 波动 幅度 较 大 。 
普通 ”私人 0.498 6.212 511.08 7 2 
普通 ”公共 3.512 11231 257.09 13 a Eero 


坦 
特殊 私人 0.642 10326 276.01 5 EE sit 
特殊 私人 0.215 8.679 735.00 7 # 2 
其 世 
蜗 


其 次 改变 病房 的 总 量 参数 进行 灵敏 度 分 析 ， 以 体现 医院 ai 
在 注重 公平 性 时 的 帕 累 托 最 优 解 的 情况 ， 改 变 医院 诊疗 单位 a 
病人 所 需 的 成 本 费用 来 表现 医院 追求 成 本 最 小 化 时 的 最 优化 2 过 由 二 
解 的 情况 。 表 7 反映 的 是 将 要 进行 的 参数 变动 。 公平 权重 指数 
表 7 参数 变动 表 图 12 三 种 不 同 算法 求 得 的 帕 累 托 最 优 解 比较 
Tab.7 Table of parameter changes Fig.12 Comparison of the Pareto optimal solutions obtained 
方案 病房 总 容量 诊疗 单位 病人 所 需 的 成 本 倍率 by three different algorithms 
280 12 在 多 目标 随机 规划 模型 中 ,设置 了 2 个 目标 函数 ， 对 应 
2 290 0.9 于 横 纵 坐标 轴 ， 分 别 为 医院 单 科室 的 运作 成 本 以 及 公平 性 权 
3 300 0.8 重 指数 。 由 于 两 个 目标 函数 之 间 存 在 关联 性 ， 在 针对 其 中 一 
4 310 0.6 个 目标 函数 进行 最 优化 时 必然 会 导致 男 一 目标 函数 的 变动 ， 
根据 上 述 的 基础 数据 设置 病人 预计 到 达 率 的 离散 度 M ”所 以 最 优 解 以 解 集 的 方式 呈现 。 


的 变化 也 会 影响 解 的 变化 ,虽然 M 较 大 时 间 题 的 近似 解 会 更 ”5.4 灵敏 度 分 析 
加 接近 于 原 多 目标 随机 规划 模型 的 求解 ， 但 是 所 耗费 的 时 间 接 下 来 分 别 对 病房 数量 以 及 医院 服务 单位 病人 所 需 成 本 
也 更 长 。 在 不 同 的 M 值 下 , 成 本 的 变动 幅度 不 大 , 但 是 公平 。 进行 调整 ， 对 病房 的 总 数量 ， 诊 疗 单位 病人 的 成 本 这 两 个 关 
权重 指数 的 变化 较 大 。 随 着 M 的 增加 , NSGA2-DS 算法 所 求 。” 键 参数 进行 灵敏 度 分 析 。 表 8 和 表 9 给 出 了 在 假设 两 点 BE! 和 
解 的 成 本 的 最 低 标 准 一 直 保 持 在 2.33 左右 ,公平 的 权重 在 0.5 ”EE? 分 别 代 表 成 本 最 小 化 和 公平 最 大 化 下 两 个 端点 的 决策 。 图 13 
左右 ， 即 在 当前 的 容量 内 可 以 接收 所 有 的 患者 。 当 M 大 于 ”和 图 14 分 别 给 出 了 三 种 算法 在 不 同 参数 变动 下 的 求解 结果 对 比 。 
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录用 定稿 陈 龙 ， 等 : 基于 改进 的 NSGA2 算法 考虑 病 患 公平 性 及 医院 运 f 究 


首先 通过 改变 病房 总 数 N 的 数量 求解 。 将 病房 总 数 设置 
为 4 个 参数 变量 280，290，300，310。 从 图 13 可 以 看 出 ， 
在 增加 N 的 同时 ， 公 平 的 比重 以 及 运作 成 本 在 不 断 地 增加 。 
病房 总 容量 的 增加 意味 着 任意 类 型 的 病人 都 可 以 获得 相对 应 
的 床位 ， 选 择 权 也 更 多 ， 同 时 容量 的 增加 也 会 提升 响应 性 以 
及 医院 病房 的 接受 率 ， 从 而 提高 公平 的 权重 ， 但 相对 应 随 着 
病床 数 的 增加 医院 的 运作 成 本 也 呈 上 升 趋势 。 总 体 上 
NSGA2-DS 求解 的 曲线 相对 于 NSGA2 与 MOBH 来 说 更 加 平 
缓 , 在 相同 的 公平 权重 指数 下 成 本 更 低 , 但 MOBH 算法 的 求 
解 所 耗 时 以 及 收敛 速度 要 略 快 于 NSGA2-DS 算法 , 在 求解 结 
果 上 不 如 NSGA2-DS， 解 集 波动 幅度 较 大 。 正 如 表 8 所 示 ， 
随 着 容量 增加 , 病房 的 接收 率 并 没有 降低 , 从 图 13 的 趋势 中 
可 以 看 出 病房 数 达 到 一 定 的 最 大 值 时 ， 两 种 算法 都 会 有 趋 近 
于 一 个 最 优 解 ， 在 这 种 情况 下 成 本 与 公平 都 会 达到 平衡 。 


237.5 


一 一 NSGA2 病 店 容量 =310 
-二 -MOBH 病房 容量 =310 
一 2 一 NSGA2-DS 病房 容量 =310 
一 * 一 NSGA2 病房 容量 =300 
琶 一 MOBH 病房 容量 =300 
232.5 一 一 NSGA2-DS 病房 容量 =300 
天- NSGA2 病房 容量 =290 
-如 -MOBH 病房 容量 =290 


一 -NSGA2-DS 病房 容量 =290 
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一 一 NSGA2 病房 容 一 =280 
一 MOBH 病房 容量 =280 


一 一 NSGA2-DS 病 店 容量 =280 
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图 13 三 种 算法 在 总 病房 数量 不 同 的 最 优 解 
Fig.13 Optimal solutions for the three algorithms with different 


numbers of total wards 
8 不同 病房 容量 下 的 最 优 帕 累 托 端点 
Tab.8 ”Optimal Pareto endpoints for different ward capacities 


病房 容量 病房 接收 率 


E! ={1,1,0.667,1,1,1,1,1,0.967,0.967,1,1,0.933,1,0.95,1,1,0.883,0.933,1} 
0 E? ={0.8,1,0.817,1,0.83,1,0.833,1,0.817,1,0.8,1,0.833,1,0.8,1,1,1,0.883,1} 
E! ={1,1,0.883,1,0.983,1,1,1,1,0.967,1,0.983,1,1,1,1,0.983,0.883,0.983,1} 
a E? ={0.85,1,1,1,0.85,1,0.883,1,0.85,1,0.867,1,0.867,1,0.867,1,0.883,1,1,1} 
El={1,1,1,1,1,1,1,L,L,1,1,1,1,1,1,1,1,1,1,1} 
0 E? ={0.9,1,1,1,0.883,1,0.883,1,0.9,1,0.917,1,0.9,1,0.917,1,1,1,1,1} 
E'={1,1,1,1,,1,1,1,1,1,L,,1,1,1,1,1,1,1,1 
310 


E? ={0.95,1,1,1,0.933,1,0.95,1,0.95,1,1,1,0.917,1,0.917,1,1,1,1,1} 


其 次 ， 在 其 他 的 参数 不 变 情况 下 ， 通 过 改变 医院 服务 各 
种 类 型 单位 病人 所 需 成 本 检验 其 对 最 优 解 集 产生 的 影响 。 将 
单位 成 本 的 倍率 依次 调整 为 1.2, 0.9, 0.8, 0.6。 如 图 14 所 示 ， 
三 种 算法 求 得 的 结果 在 公平 权重 指数 波动 时 运作 成 本 的 变化 
幅度 明显 小 于 改变 病房 容量 下 所 求 结 果 ， 大 多 数 情况 下 


所 求 结果 最 优 ， 其 次 ， 在 其 他 的 参数 不 变 情况 下 ， 通 过 改变 
医院 服务 各 种 类 型 单位 病人 所 需 成 本 检验 其 对 最 优 解 集 产生 
的 影响 。 将 单位 成 本 的 倍率 依次 调整 为 1.2，0.9，0.8，0.6。 
如 图 14 所 示 ， 三 种 算法 求 得 的 结果 在 公平 权重 指数 波动 时 
运作 成 本 的 变化 幅度 明显 小 于 改变 病房 容量 下 所 求 结果 ， 大 
多 


数 情 况 下 NSGA2-DS 所 求 得 成 本 在 相同 公平 指数 下 低 于 
他 两 种 算法 ， 所 求 结果 最 优 ， 解 的 稳定 性 较 强 ， 但 在 收敛 
速度 上 略 慢 于 MOBH 算法 。 从 表 9 和 图 14 可 以 看 出 在 所 有 
最 优 解 中 E'! 和 至? 的 决策 是 相同 的 ， 在 成 本 呈现 倍数 变化 
同时 ， 决 策 变量 略微 受到 影响 且 不 如 病房 总 量变 化 对 求解 


NSGA2-DS 所 求 得 成 本 在 相同 公平 指数 下 低 于 其 他 两 种 算法 ， 


结果 的 
增加 而 减少 。 由 于 运作 成 本 是 由 服务 单位 病人 所 需 成 本 与 
人 和 总数 相 乘 ， 同 时 病人 总 数 与 医院 的 接受 率 有 关 ， 因 此 总 
成 本 会 随 着 单位 成 本 的 减 小 而 降低 , 而 公平 权重 指数 基 
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的 影响 大 ， 帕 累 托 边界 上 的 最 优 解 的 数量 随 着 成 本 的 


|e 


斌 


保持 不 变 。 


表 9 单位 成 本 不 同 倍率 下 的 最 优 帕 累 托 端点 


Tab.9 Optimal Pareto endpoints at different multipliers of unit cost 


倍率 病房 接收 率 

0.6 盏 = 人 LLLLLLLLLLLLLLLLLLIT 

0.8 盏 = 人 LLLLLLLLLLLLLLLLLLIT 

0.9 EE?={0.9,1,1,1,0.88333,1,0.88333,1,0.9,1,0.91666,1,0.9,1,0.91666,1,1,1,1,1} 
1.2 E?={0.9,1,1,1,0.88333,1,0.88333,1,0.9,1,0.91666,1,0.9,1,0.91666,1,1,1,1,1} 


本 文 在 对 多 目标 随机 规划 模型 进行 线性 化 成 混合 整数 线 


性 规划 的 基础 上 采用 NSGA2-DS 算法 进行 求解 ， 为 了 验证 
NSGA-DS 算法 的 有 效 性 以 及 前 沿 性 ,与 原 有 的 NSGA2 算法 ， 


以 及 引入 的 MOBH 算 法 进行 比较 ,图 13 和 图 14 显示 NSGA2- 
DS 算法 求 得 解 都 优 于 其 他 两 种 算法 。 三 种 算法 所 得 求 最 优 


解 的 解 距 随 着 病房 总 数 的 增加 而 减少 。 根 据 图 13 与 图 14 的 
趋势 分 析 ， 当 病房 数量 达到 一 定 的 最 大 值 时 候 ， 三 种 算法 都 
会 求 得 最 优 解 ， 使 得 成 本 与 公平 达到 平衡 ， 而 单位 成 本 对 于 


整体 运作 成 本 的 影响 较 大 ， 对 公平 性 几乎 没有 影响 。 


6 


为 影响 
喘 反映 


性 最 佳 
划 模 型 
便 求 解 
复杂 问题 上 全 局 搜索 能 力 较 差 、 收 敛 性 不 高 、 分 布 性 较 差 的 
局 限 性 进行 改进 ， 设 计 改 进 后 的 NSGA2-DS 算法 。 最 后 先 通 


300 


-AA 一 万。 | 一人- NSGA2 倍率 =1.2 

全 一 和 一 MOBH 信 素 =1.2 
A 一 各 NSGA2-DS 倍率 =1.2 

mg 2 本 一 #_NSGA2 倍率 -0.9 

一 一 MOBH 信 率 -0.9 
一 9 一 NSGA2-D5 倍率 -0.9 

260 一 一 NSGA2 信和 率 -=0.8 

一 一 MOBH 倍率 =0.8 


站 一 一 N5GA2-D5 倍率 =0.8 
EE 240 SD 一 e 一 NSGA2 倍率-0.6 
长 一 人 4 一 9 一 MOBH 售 率 =0.6 
: 一 一 %— NS5GA2-D5 倍率 =0.6 
惧 


公平 权重 指数 


图 14 三 种 算法 在 单位 成 本 倍率 不 同 的 最 优 解 


Fig.14 Optimal solutions for three algorithms 


with different unit cost multipliers 


结束 语 


本 


全 到 


文 提出 了 一 种 以 每 位 患者 所 能 接受 的 最 低 等 待 天 数 作 
姑 素 ， 考 虑 准 入 性 与 响应 性 公平 的 权重 测度 指标 来 反 
医院 与 患者 的 公平 性 ， 并 结合 医院 运作 成 本 以 及 公平 
为 目标 建立 多 目标 随机 规划 模型 。 针 对 多 目标 随机 规 
的 复杂 性 ， 采 用 线性 化 将 原 有 复杂 模型 转换 转换 成 方 
的 混合 整数 线性 化 模型 。 对 NSGA2 算法 在 求解 高 维 


过 实验 检验 了 NSGA2-DS 算法 的 性 能 , 算法 在 测试 函数 和 文 


章 的 问 
不 同 病 
度 分 析 
模型 的 有 效 性 以 及 适用 性 。 


本 


题 方面 效果 均 优 于 NSGA2 与 MOBH 算法 ， 其 次 调整 
房 容量 、 服 务 不 同 病 人 下 的 单位 成 本 等 参数 进行 灵敏 
， 得 出 参数 的 变化 与 目标 之 间 的 关联 ， 进 一 步 验 证 了 


文 研究 的 方向 主要 是 针对 单 科 室 的 病床 分 配 ， 未 来 的 


研究 方向 可 以 进一步 的 拓展 ， 一 方面 考虑 完全 公平 情况 下 ,上 
考虑 私人 病房 预 留 一 定 的 床位 给 VIP 病人 ， 同 时 病人 的 类 型 更 


“ey 


ul 


丰富 ,病人 的 行为 和 等 待 时 间 的 复杂 关系 等 等 , 这 是 难点 所 在 ; 


另 一 方面 不 仅仅 局 限于 单 科 室 ， 可 以 考虑 复杂 多 科室 乃至 医院 


集群 


亡 


面 的 临床 资源 联合 优化 ， 为 实际 的 医疗 资源 提供 更 多 参考 。 
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